arXiv:1506.05770vl [math.OC] 18Jun2015 


Distributed Verification of Structural Controllability for 

Linear Time-Invariant Systems 

Joao Carvalho Sergio Pequito A. Pedro Aguiar 
Soummya Kar ^ Karl H. Johansson § 


Abstract 

Motivated by the development and deployment of large-scale dynamical systems, often composed of geo¬ 
graphically distributed smaller subsystems, we address the problem of verifying their controllability in a distributed 
manner. In this work we study controllability in the structural system theoretic sense, structural controllability. In 
other words, instead of focusing on a specific numerical system realization, we provide guarantees for equivalence 
classes of linear time-invariant systems on the basis of their structural sparsity patterns, i.e., location of zero/nonzero 
entries in the plant matrices. To this end, we hrst propose several necessary and/or sufficient conditions to ensure 
structural controllability of the overall system, on the basis of the structural patterns of the subsystems and their 
interconnections. The proposed verification criteria are shown to be efficiently implementable (i.e., with polynomial 
time complexity in the number of the state variables and inputs) in two important subclasses of interconnected 
dynamical systems: similar (i.e., every subsystem has the same structure), and serial (i.e., every subsystem outputs 
to at most one other subsystem). Secondly, we provide a distributed algorithm to verify structural controllability 
for interconnected dynamical systems. The proposed distributed algorithm is efficient and implementable at the 
subsystem level; the algorithm is iterative, based on communication among (physically) interconnected subsystems, 
and requires only local model and interconnection knowledge at each subsystem. 

1. Introduction 

Recent years have witnessed a growing use of large-scale dynamical systems, notably, those with a 
modular structure [[U, [O, such as content delivery networks, social networks, robot swarms, and 

smart grids. Such systems, often geographically distributed, are composed of smaller subsystems, which 
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we may refer to as agents. A typieal eoneern is that of ensuring that the system, as a whole, performs as 
intended. More than often, while analyzing these intereonnected dynamical systems, which in this paper 
are considered as continuous linear-time invariant (LTI) subsystems, the exact parameters of the plant 
matrices are unknown. Therefore, we focus on the zero/nonzero pattern of the system’s plant, which we 
refer to as sparsity pattern, and, in particular, on structural counterpart of controllability, i.e., structural 
controllability Hi. 

It is worthwhile noting that these agents may be homogeneous or heterogeneous, from the structural 
point of view. When the agents are homogeneous, their plants and connections have the same sparsity 
pattern and the system is referred to as a similar system. Otherwise, the agents are heterogeneous and two 
possible scenarios are conceivable: (i) an agent may receive input from (possibly several) other agents 
but it only feeds to one other agent, in which case the overall system is referred to as serial, and {ii) the 
interactions between agents can be arbitrary, typically expected in setups where the inputs to the agents is 
broadcasted from the others. In this paper, all the above subclasses of interconnected dynamical systems 
are of interest and explored in detail. More precisely, several necessary and/or sufficient conditions are 
proposed to ensure key properties of the system; further, these properties can be verified by resorting to 
efficient algorithms, i.e., with polynomial time complexity in the number of state and input variables. 

In some applications, the problem of composability is particularly relevant; for example, a swarm of 
robots possessing similar structure where the interaction (possibly through communication or information 
exchange) topology may change over time, or where robots may join or leave the swarm over time. Then, 
the existence of necessary and/or sufficient conditions on the structure and interconnection between these 
agents contribute to controllability-by-design schemes, i.e., we ensure that by inserting an agent into the 
interconnected dynamical system the dynamical system continues to be controllable. In addition, we can 
also specify with which agents an agent should interact with such that those conditions hold. 

A swarm of robots can also be composed of a variety of heterogeneous agents in which case controllability- 
by-design is also important. Yet, due to constraints on the communication range, the interaction between 
agents is merely local, even if some additional information is known. Therefore, in the context of serial 
systems we can equip each subsystem with the capability of inferring if the entire system is structurally 
controllable or not. In other words, distributed algorithms that rely only on interactions between a 
subsystem and its neighbors, where information about their structure may be shared. In particular, if we 
equip the robots in the swarm with actuation capabilities that can be activated when the interconnected 

*The sparsity pattern (or structural matrix) of a real matrix M is a binary matrix M, satisfying Mij = 0 if and only if Mi^j = 0. We 
then call the pair of structural matrices {A, B) associated to an LTI system, a structural system, and say that it is structurally controllable 
if and only if there exists a controllable pair of real matrices {A',B') with the same sparsity pattern as {A,B). In fact, we notice that if 
{A, B) is structurally controllable, then almost all pairs of real matrices (,4, B) with the same sparsity pattern as {A, B) are controllable |4(. 
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dynamical system is not structurally controllable, we can render this interconnected dynamical system 
structurally controllable. 

Nonetheless, imposing a priori knowledge of the structure of the interconnections in the system (for 
instance, whether it is a serial system) can be restrictive, so distributed algorithms that aim to verify 
structural controllability of general interconnected dynamical systems are desirable. Hereafter, we provide 
such an algorithm: it requires the interaction between a subsystem and its neighbors, but it does not require 
(directly) sharing the structure of the subsystems involved. Instead, it requires only partial information 
about its structure, which is also desirable from a privacy viewpoint. The proposed scheme is also 
particularly suitable to other applications such as the smart grid of the future, that consists of entities 
described by subsystems deployed over large distances; in particular, notice that in these cases, the different 
entities may not be willing to share information about their structure due to security or privacy reasons. 

Related Work: Structural controllability was introduced by Lin [j5l in the context of single-input 
single-output (SISO) systems, and extended to multi-input multi-output (MIMO) systems by Shields and 
Pearson 0. A recent survey of the results in structural systems theory, where several necessary and 
sufficient conditions are presented, can be found in Q. 

In this paper, we are interested in understanding how the connections between different dynamical 
subsystems enable or jeopardize the structural controllability of the overall system. This study follows the 
general lines provided in [17]|, but the verification procedures proposed were based in matrix nets leading 
to a computational effort that increases exponentially with the dimension of the problem. Alternatively, 
an efficient method is proposed in [[3l that accounts for the whole system instead of local properties (i.e., 
the components of the system and their interconnections), but this method does not apply to arbitrary 
systems. More precisely, it is assumed that when connected, the state space digraph (see Section 
for formal definition) is spanned by a disjoint union of cycles, which is called a rank constraint. In 
contrast, in flU and [(9l, the authors have presented results on the structural controllability of interconnected 
dynamical systems with interconnection in cascade. Nevertheless, these structures are not unique, and 
the interconnection of these is established assuming such connectible structures are given; therefore, no 
practical criteria for computing such structures and verifying the results were provided. More recently, 
in [[Toll similar results were obtained by exploring which variables may belong to a structure and referred 
to as controllable state variables. Thus, similarly to IH and [|9||, the results depend on the identified 
structures, but no method to systematically identify these structures is provided. In [fTTll the study is 
conducted assuming that all the subsystems, except a central subsystem which is allowed to interact with 
every other subsystem, have the same dynamic structure and the interconnection between the several 
subsystems also has the same structure (even though they may not be used). The implications of local 
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interactions on the system’s structural controllability determined in ifTTll can be obtained as particular 
solutions to the results proposed hereafter. 

In [[T2|, we studied the problem of determining the sparsest input matrix ensuring structural controlla¬ 
bility, and additionally addressed the sparsest co-design problem, i.e., determining the sparsest solution of 
inputs, outputs and information patterns such that decentralized control laws guaranteeing that the closed- 
loop system enables arbitrary pole placement, based in static output feedback, exist. Further, polynomial 
algorithms with computational complexity 0{n^) were provided to both problems, where n is the the 
number of state variables. In [fT3ll . we extended the results in [|T^ . to the setting where the selection 
of inputs is constrained to a given collection of inputs, and sparsity of the inputs is not accounted for. 
And later, in llT4l the problem in ifT^ was further extended to determining the input matrix incurring 
in the minimum cost with respect to some cost function, and ensuring structural controllability. Further, 
procedures with 0{n‘^) computational complexity where provided, where u < 2.373 is the lowest known 
exponent associated with the complexity of multiplying two n x n matrices. All these contrast with the 
one addressed in the current paper in the sense that we aim to verify structural controllability properties 
in a distributed fashion. 

On the other hand, composability aspects regarding controllability have been heavily studied by several 
authors, see for instance, [fTSlI . [fT^ . [fT71 . [fTSl . l(T9ll . [|20]|. Briefly, all these resort to the well known 
Popov-Belevotch-Hautus (PBH) eigenvalue controllability criterion for LTI systems ll2T]| . We notice that 
this criterion requires the knowledge of the overall system to infer its controllability. The reason is closely 
related with the loss of degrees of freedom imposed by interconnected dynamical systems, as well as 
conservation laws in general, that result in the reduction of the rank of the system’s dynamics matrix 
when compared with the sum of the ranks of the dynamics matrices of each subsystem. Consequently, 
even if all subsystems are controllable, after their interconnection the resulting dynamical system may 
not be. Notwithstanding, the same does not happen when dealing with structural systems, where if all 
subsystems are structurally controllable, then the overall system is structurally controllable. This is why, the 
key technical results in this paper correspond to verifying structural controllability in scenarios where the 
individual subsystems are not structurally controllable by themselves. In addition, while not guaranteeing 
that a system is controllable, we can regard the structural controllability analyses in this paper as necessary 
conditions for controllability. 

Main Contributions: The main contributions of this paper are threefold: 

(i) We provide sufficient conditions for similar systems to be structurally controllable. More precisely, 
these rely only on the structure of the subsystem and interconnection between subsystems. A distributed 
algorithm is proposed, that can verify these conditions in polynomial time. 
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(ii) We provide sufficient conditions for serial systems to be structurally controllable. A distributed 
algorithm to verify these conditions is provided. It requires only the knowledge of the subsystem and 
its neighbors’ structure, as well as its interconnections. This algorithm requires that a subsystem has the 
capability to interact with its neighbors only, and has computational complexity equal to O (nl), where 
ns denotes the total number of state variables and inputs present in a subsystem and its neighbors. 

(in) We provide a distributed algorithm to verify necessary and sufficient conditions to ensure structural 
controllability for any interconnected dynamical system that consists of LTI subsystems. This algorithm 
requires that a subsystem has the capability to interact with its neighbors only, have access to its own 
structure and partial information regarding decisions performed by its neighbors that do not require sharing 
(directly) the structure of the neighboring agents. 

The rest of this paper is organized as follows. In Section we formally describe the problem statement. 


Section III introduces some concepts in structural systems theory, that will be used throughout the 
remainder of the paper. The main contributions are presented in Section IV and in Section |V] we provide 
examples that illustrate the main findings. Finally, Section |V^ concludes the paper and discusses avenues 
for further research. 


II. Problem Statement 

Consider r linear time-invariant (LTI) dynamical systems described by 

Xi(t) = AiXiit) + BiUiit), i = 


where Xi G M”* is the state, and Ui G the input of the i-\h system. The dynamical system can be 
described by the pair (Aj, Bi), where Aj G is the dynamic matrix of subsystem i and Bi G 

its input matrix. By considering the interconnection from subsystem i to subsystem j for all possible 
subsystem pairs we obtain the interconnected dynamical system described as follows: 
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where the state is given by x = ... g with n = Y7i=i the input given by m = 

[ul ... ul]^ G with p = addition, Ei^ G is referred to as the connection matrix 

from the j-th subsystem to the i-th subsystem. Further, we denote the system ([T]) by the matrix pair 
(A, B), denoting the i-th subsystem, i = 1,..., r of ([T]) by the matrix pair (Aj, Bi). Finally, we call those 
subsystems (Aj, Bj), with j = 1,..., r such that Ej^i ^ 0, the outgoing neighbors of the i-th subsystem. 
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and those that Eij ^ 0 the incoming neighbors of the i-th subsystem; we refer to them eolleetively as 
the neighbors of the i-th subsystem. 

Now, eonsider the sparsity pattern of the matrix pair {A, B) whieh we denote by the struetural system 
(A, 5); similarly, we denote by {Ai, Bi) the struetural pair of matriees assoeiated with {Ai,Bi), and 
Eij the sparsity pattern of Ej^. Within the eontext of struetural systems, the struetural eounterpart of 
eontrollability ean be introdueed as follows. 

Definition 1 (^): Given a struetural system {A,B), we say that it is strueturally eontrollable if and 
only if, there exists at least one eontrol system {A, B) with the same sparsity pattern as {A, B) (i.e. Aij = 0 
if Aij = 0 and Bi^k = 0 if Bi ^ = 0) whieh is eontrollable. 

It ean be seen, from density arguments, that if (A, B) is strueturally eontrollable, then almost all 
eontrol systems {A^B) with the same sparsity as {A,B) are strueturally eontrollable H. We say that a 
eontrol system {A, B) is structurally controllable if the assoeiated struetural system {A, B) is strueturally 
eontrollable. 

Therefore the problem addressed in the eurrent paper ean be posed as follows. 

Problem: Given a eolleetion of eontrol systems (Ai,Bi), i = 1,... ,r, and the intereonneetion from 
eaeh subsystem i to its neighbors, i.e., {Aj, Bj, Ejj) for all j ^ i, design a distributed proeedure to 
determine if the intereonneeted eontrol system (A, B) given in ([1]) is strueturally eontrollable. o 

Further, note that in a non-struetural setting, loeal properties are not enough to guarantee eontrollability 
sinee the eonneetion to other subsystems may lead to parameter eaneellations ffT9l . [|20l : therefore, by 
resorting to struetural system theoretie study, the approaeh presented hereafter allows us to obtain only 
neeessary eonditions for eontrollability. 

III. Preliminaries and Terminology 

In this seetion, we review some preliminary eoneepts used to analyze the problem of struetural eon¬ 
trollability of intereonneeted dynamieal systems. 

A. Structural Systems 

In order to perform struetural analysis effieiently, it is eustomary to assoeiate to ([1]) a direeted graph, 
or digraph V = {V,S), in whieh V denotes the set of vertices and S the set of edges, where {vj,Vi) 
represents an edge from the vertex Vj to the vertex Vi. To this end, let A G {0,and B G {0, 
be the binary matriees that represent the sparsity patterns of A and B as in ([T]), respeetively. Denote 
by A" = {xi, ■ ■ ■ ,Xn} and U = [ui, ■ ■ ■ ,Up} the sets of state and input vertices, respeetively. And by 
= {{xi,Xj) : Aji 7 ^ 0}, = {{uj,Xi) : Bij 7 ^ 0}, the sets of edges between the vertex 

sets in subseript. We may then introduee the state digraph 'D{A) = {?C,Ex,x) and the system digraph 


7 


V{A,B) = {X UU,Sx,x £u,x)- Note that in the digraph V{A,B), the input vertices representing 
the zero columns of B correspond to isolated vertices. As such, the number of effective inputs, i.e., the 
inputs which actually exert control, is equal to the number of nonzero columns of B, or, in the digraph 
representation, the number of input vertices that are connected to at least one state vertex through an edge 
in ^u,x- 

A directed path from the vertex vi to is a sequence of edges {(^ 1 ,^ 2 ), {^ 2 X 3 )^ ■ ■ ■ i {vk-i,Vk)}. If 
all the vertices in a directed path are distinct, then the path is said to be an elementary path. A cycle is 
an elementary path from vi to Vk, together with an edge from Vk to Vi. 

Given a digraph V = (V, S), we say that V' = (V', S') is a subgraph of V if it is a digraph with V' C V 
and S' C S, which we denote by T)' C T). Further, we say that T)' spans "D if = V. 

We also require the following graph-theoretic notions [|22l : A digraph V is strongly connected if there 
exists a directed path between any two vertices. A strongly connected component (SCC) is a subgraph 
= (Vs, Ss) of V such that for every m, n G Vs there exist paths from m to f and from n to m and is 
maximal with this property (i.e., any subgraph of V that strictly contains Vs is not strongly connected). 

Definition 2 t/l72l/).' An SCC is said to be linked if it has at least one incoming or outgoing edge from 
another SCC. In particular, an SCC is non-top linked if it has no incoming edges to its vertices from the 
vertices of another SCC. o 

Further, given a digraph V = {V,S), we say that "D is a weakly connected digraph if {V,S U S'^) is 
strongly connected, where S'^ = {(V,n) : {v,v') e S} 

For any digraph V = (V, S) and any two sets 5i, ^2 C V we define the bipartite graph B{Si, S 2 , ^^ 51 , 53 ) 
where we call the set of left vertices, and S 2 the set of right vertices', and the edge set Ss^^s 2 = 

n (iSi X ^ 2 ). We call the bipartite graph B{V, V, S) the bipartite graph associated with V{V, S). In the 
sequel we will make use of the state bipartite graph, B{A) = B{X, X, Sx,x), which is the bipartite graph 
associated associated with the state digraph V{A) = {X, Sx,x), and the system bipartite graph B{A, B) = 
B{U U X,X,Sx,x U Su,x)- An illustration of the state bipartite graph can be found in Figure [^(&) 
associated with the state digraph in Figure (a). 

Given a bipartite graph B{S\,S 2 ,Ss^^s 2 )^ a matching M corresponds to a subset of edges in Ssi,s 2 
so that no two edges have a vertex in common, i.e., given edges e = (si,S 2 ) and d = (s'i,S 2 ) with 
si, s'l e and S 2 , s '2 G S 2 , e, e' G M only if si s'l and §2 s' 2 . An example of a matching is provided 

in Figure [^(6). 

A maximum matching M* is a matching M that has the largest number of edges among all possible 
matchings, as the one shown in Figure [^(6). 

Further, it is possible to assign a weight to the edges in a bipartite graph, say c(e) (where c is a function 


(a) (b) M* 



Fig. 1: In (a) we have depicted a digraph 'D{A). The SCCs are depicted within the dashed boxes; the digraph consists of two non-top 
linked SCCs ,^ 2 ) and two non-bottom linked SCCs (A/j®, A/”^). In (6) we depicted the state bipartite graph B{A) = B{X, X, £x,x), 
along with a maximum matching M* comprising the edges in red, which are also highlighted in (a) as they comprise a minimum path and 
cycle decomposition, as explained in Lemma 

from £ 5^,82 to We thus obtain a weighted bipartite graph, and can introduce the concept of minimum 
weight maximum assignment problem. This problem consists in determining a maximum matching whose 
overall weight is as small as possible, i.e., a matching M'^ such that 

= arg min > c(e), 

M&M ^ 

eSM 

where M. is the set of all maximum matchings. This problem can be efficiently solved using the Hungarian 
algorithm ^23t . with complexity of O (max{|5i|, |52|}^). We call the vertices in iSi and S 2 belonging to 
an edge in M*, the matched vertices with respect to (w.r.t.) M*, otherwise, we call them unmatched 
vertices. It is worth noticing that there may exist more than one maximum matching, for example, in 
Figure [^(6) replacing the edge {x 8 ,xg) in M* with {xg,xg) yields a different matching with the same 
number of edges. For ease of referencing, in the sequel, the term right-unmatched vertices, with respect 
to 1 S 2 ,£^ 51 , 52 ) a matching M, not necessarily maximum, will refer to those vertices in S 2 that 

do not belong to an edge in M*, and by duality a vertex from iSi that does not belong to an edge in M* 
is called a left-unmatched vertex. 

The following result gives us some liberty in choosing the right- and left-unmatched vertices of the 
maximum matching of a bipartite graph. 

Lemma 1 1/1721/ ).' Let B{V,V,£) be a bipartite graph. If and are two possible maximum 
matchings of B{V,V,£) with right-unmatched and left-unmatched vertices given by {Ur{M^),Ul{M^)) 
and {Ur{M‘^),Ul{M‘^)) respectively, then, there exists a maximum matching M* of B{y,V,£) with sets 
of right-unmatched and left-unmatched vertices given by {Ur{M^), Ul{M‘^)). o 

Now, we can interpret a maximum matching of a bipartite graph associated to a digraph, at the level 
of the digraph. 

Lemma 2 (Maximum Matching Decomposition /f72]/ ).- Consider the digraph V = {y,£) and let M* be 
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a maximum matching associated with the bipartite graph Then, the digraph V = (V,M*) 

comprises a disjoint union of cycles and elementary paths (by definition an isolated vertex is regarded 
as an elementary path with no edges), beginning in the right-unmatched vertices and ending in the left- 
unmatched vertices of M*, that span V. Moreover, such a decomposition is minimal, in the sense that no 
other spanning subgraph decomposition of T>{A) into elementary paths and cycles contains strictly fewer 
elementary paths. o 

In addition, to make comparisons with previous work (namely, O and (iH), we need to introduce the 
following definitions. 

Definition 3 f/E]/)-' Given a digraph V, an elementary path in V, also called a stem, is a cactus. Given a 
cactus Q = {Vg,Sg) C V, and a cycle C = {Vc,£c) ^ such that Q and C have no vertices in common, 
and there is an edge from a vertex in ^ to a vertex in C, then Q UC = {Vg U Vc, Eg U Eg) is a cactus, o 

Particularly, in the case where V = V{A,B), a cactus ^ in P is called an input cactus if the stem 
starts at an input vertex. Further, we note that the decomposition into disjoint elementary paths and cycles 
in Lemma 1^ can be used to determine a spanning decomposition of the graph into disjoint cacti [|T^ . 

When dealing with interconnected dynamical systems, the structure of the connection between the 
subsystems will create connections between the SCCs of different subsystem digraphs. This, in turn, 
makes it difficult to identify the SCCs of the system digraph of the overall system by analysing the 
SCCs of each subsystem digraph seperately and the connections to their neighbors. Hence, we introduce 
the concept of reachability [31. We say that a state vertex x in a system digraph is input-reachable or 
input-reached if there exists a path from an input vertex to it. 

All of these constructions can be used to verify the structural controllability of an LTI system by 
analysing the associated graphs. 

Theorem 1 (/3]/, For LTI systems described by ([T]), the following statements are equivalent: 

(1) The corresponding structured linear system {A, B) is structurally controllable; 

(2) The digraph V{A,B) is spanned by a disjoint union of input cacti; 

(3i) The non-top linked SCCs of the system digraph T>{A, B) are comprised of input vertices, and 
(3 m) there is a matching of the system bipartite graph B{A, B) without right-unmatched vertices; 

(4i) Every state vertex is input-reachable, and 

(4m) there is a matching of the system bipartite graph B{A, B) without right-unmatched vertices. 

o 

IV. Main Results 

We begin this section by providing sufficient conditions for an interconnected dynamical system to be 
structurally controllable in the case where all the subsystems have the same structure (Theorem and 
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Theorem [^. We then foeus on more general interconneeted dynamieal systems, ealled serial systems, 
and propose suffieient eonditions for their structural controllability (Lemma [^; in addition, an efficient 
distributed algorithm (Algorithm to verify these conditions is provided, which has its correctness and 
complexity proven in Theorem In light of these conditions, we explain why previous results in this 
line dSl presented conditions that are only sufficient instead of necessary and sufficient (Figure [^. Finally, 
we end this section by describing an efficient distributed algorithm (Algorithm to verify the structural 
controllability of an arbitrary interconnected dynamical system, which has its correctness and complexity 
proven in Theorem In order to perform this verification, each subsystem has to perform calculations 
using the information about itself and its neighbors. 

Often, interconnected dynamical systems under analysis are composed of subsystems that are similar, 
as formally introduced next. 

Definition 4: Let E G {0,1}’’^^, A',H G {0,1}”^”, B' G {0,1}”^^, be matrices with the restriction 
that Ei^i = 0 {i = l,...,r). Then, we denote by {A',B',H,E') the structural system (A,B), where 
A = A') y {E ® H) and B = Ir® B', where ® denotes the matrix Kronecker product. Further, we 

say this system to be composed of r similar subsystems, or, in short, a similar system. o 

Remark 1: Note that in the case of similar systems, H is the structural matrix modeling the interactions 
between each subsystem and its neighbors, all of which have the same structure. 

Definition 5: Let {A, B) be the structural matrices associated with the interconnected dynamical system 
in O- We define the condensed graph of the system as being the digraph V*{A) = V{A,S), where 
Oi e A = {oi,..., Or} is a vertex representing the i-th subsystem, and (a*, af E S = {(a*, aj)\Ej^i A 0} 
a directed edge representing a communication from subsystem j to subsystem i. Moreover, if there is no 
directed edge ending in a vertex, this vertex is referred to as a source. o 

Note that in the case that a system {A, B) is composed of r similar systems and parametrized by 
matrices {A', B', H, E) the condensed graph V*{A) is the same as the digraph D{E). Now, we assess the 
structural controllability of these systems when their subsystems are not structurally controllable. 

Theorem 2: Let the system {A, B) be composed of r similar components, and parametrized by (A', B', H', E), 
where {A',B') is not structurally controllable, and B{A',B) has a matching without right-unmatched 
vertices. The pair (A, B) is structurally controllable if and only if (A' V H, B') is structurally controllable 
and the condensed graph V*{A) has no sources. o 

Proof: To prove the equivalence, we begin by proving that the conditions are sufficient by contra¬ 
positive; subsequently, we prove directly that the conditions are also necessary. 

First, notice that since (A', B') is not structurally controllable despite B{A', B') having a maximum 
matching without right-unmatched vertices, it follows, from Theorem [^(4), that V{A',B') has a vertex 
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which is not reachable from any input vertex. Subsequently, assume that V*{A) has a souree, then it 
follows that there is a subsystem (^4' , 5') with no ineoming edges from other subsystems, and so the 
overall system digraph V{A^B) has a state vertex without a path from any input vertex to it. Henee, by 
Theorem[^(4), {A, B) is not strueturally eontrollable. Further, {A'\/H, B') is not strueturally eontrollable. 
Sinee B{A',B') has a matehing without right-unmatehed vertiees, so does B{A' V H,B'), whieh means, 
by Theorem [^(4) that, V{A' V H, B') must have a state vertex whieh is not reaehable from any input 
vertex. This implies that the eorresponding state vertex is not reaehable from an input vertex in any of 
the subsystems (sinee a path from an input vertex in the overall system translates into one sueh path in 
V{A'y H,B')). 

Finally, assume that {A' V H,B') is strueturally eontrollable and that V*{A) has no sourees. Then, for 
eaeh state vertex, there is a path from an input vertex to it, whieh implies, by Theorem [^(4), that {A, B) 
is strueturally eontrollable. ■ 

In the next result, we relax the assumptions from Theorem about the strueture of the dynamies of 
the subsystems. Thus, allowing for applieations in the design of intereonneetions between subsystems that 
may fail to meet sueh eriteria. 

Theorem 3: Given an intereonneeted dynamieal system (^4, B) eomposed of r similar eomponents, and 
parametrized by (A', B', H, E), where {A',B') is not strueturally eontrollable, then (A, B) is strueturally 
eontrollable if both {A' V H,B') is strueturally eontrollable and I)*(A) is spanned by eyeles. o 

Proof: First, notiee that if the digraph V*{A) is spanned by eyeles, every vertex belongs to a eyele, 
and, in partieular, it means that T)*{A) has no sourees. In this ease, the method of proof of Theorem]^ 
is applieable to show that every state vertex has a path from an input vertex to it, so all that remains to 
show is that the B{A,B) has no right-unmatehed state vertiees with respeet to some maximum matehing. 
To this end, we first assume (without loss of generality) that I)*(A) has one spanning eyele, and that the 
subsystems {Af B[), • • • , [Af B'A are ordered in sueh a way that = 1 for i = 1,... ,r — 1, and 

El^r = 1 - 

Now, denote the state and input vertiees of the z-th subsystem by x\ with k = I,-- - ,?t, and u\ 
with I = 1, • • • ,m, respeetively. In addition, let M be a maximum matehing of B{A' V H,B') without 
right-unmatehed state vertiees, then we ean partition M' into three matehings M'^, Mf, M'^ eomprising, 
respeetively. The edges in M' are of the form {ui,Xk), {xi,Xk) when Ak^i = 1, and the remaining ones 
are of the form {xi,Xk), when A^^i = 0 and = 1. Finally, eonsider the matehing M of B{A,B) 
eomprising the edges: 

. iui,xi), if {Uk,xi) is in 
. ixi,xi), if {xk,xi) is in M'^-, 
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. {xl,x\), if {Xj,xi) is in M'^-, 

• {xl,x\^^), if {Xj,xi) is in M'^. 

To show that M has no right-unmatched vertices associated with it, consider a state vertex x\. of B), 
and since M' has no right-unmatched vertices, Xk is not right-unmatched in thus, there is an 

edge {xi, Xk) for some I or {ui'^Xk) in M', but this implies that either {xlxi), {x\ or (4,4) e M 

(where i — l = r when i = 1) from the construction above, thus B{A^ B) has no right-unmatched vertices 

w.r.t. M. 


Finally, if more than one cycle is necessary to span the digraph V*{A), then the same argument applies 
to each cycle individually. ■ 

Remark 2: The conditions in Theorem]^ are not necessary, consider the example in Figure]^ (6), that 
is structurally controllable, yet the condensed graph is not spanned by cycles. 


ia)V{A',B') ic)V*{A,B) 



(b) V{A,B) M 



Fig. 2: In (a) we provide the digraph T>{A',B') of a system that is not structurally controllable, where we represent in blue the single 
input vertex. By connecting three of these systems together as in (6) the system {A, B) becomes structurally controllable, as evidenced by 
the fact the matching M, associated to the path and cycle decomposition (Lemma [TJ depicted by the red edges, has no right-unmatched 
state vertices (since every state vertex has an incoming red edge), and the fact that every non-top linked SCC is comprised of an input 
vertex. Finally, in (c) we show that the condensed graph B*{A) of the system in (6) is not spanned by cycles, showing that the condition 
in Theorem is not necessary. 

Now, we move toward methods for verifying struetural eontrollability of intereonneeted dynamieal 
systems by resorting to distributed algorithms. To this end, we begin by introdueing a result that will 
allow us to infer struetural eontrollability of a family of intereonneeted dynamieal systems that we refer 
to as serial systems. Later, we provide a eomputational method to perform this verifieation in a distributed 
manner, that is, in whieh eaeh subsystem only needs to have partial information about the system in order 
to verify if the system is strueturally eontrollable or not. 

It is worth noting that in order to ensure that all of the algorithms in the present paper ean work as 
intended, several assumptions have to be made about the subsystems, and how they internet. Namely, that 
eaeh subsystem has a proeessing unit, and ean send arbitrary messages to its neighboring subsystems; in 
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addition, each subsystem is aware of the number of subsystems in the overall system and possesses a 
unique id, and that the condensed graph of the system V*{A) is weakly connected. 

Lemma 3: Consider the structural system {A, B) as in o with subsystems {Ai, Bi ),..., {A^, B^). 
Then the system {A, B) is structurally controllable, if there exist maximum matchings Mq, ..., Mr of the 
bipartite graphs B{Ai ),..., B{Ar) respectively, such that the following conditions hold: 

1) For each subsystem {Aj, Bj), j = 1,..., r, the non-top linked SCCs of 'D{Aj, Bj) consist of input 
vertices; and, 

2) the bipartite graph 



admits a perfect matching, where UiiMi) and UR^Mi) are the sets of left- and right-unmatched 
vertices of Mi respectively, and £uL(Mj),UR{Mi) A Bx,x is the set of edges from vertices in UL{Mj) 


to vertices in UR^Mi). 


o 


Proof: First, note that the non-top linked SCCs of V{A, B) consist of SCCs of the subsystem digraphs 
V{Ai,Bi). Further, for one such SCC to be non-top linked, it must contain at least one non-top linked 
see of some V{Ai, Bi) in it. Therefore, since every non-top linked SCC of every V{Ai, Bi) is comprised 
of input vertices, and there are no edges from any neighboring system to input vertices, the non-top linked 
SCCs of 'D{A, B) must contain input vertices. 

Secondly, note that the union of the maximum matchings Mi of the B{Ai, Bi) comprises a matching M 
of B{A,B). Further, let M' be the matching mentioned in the second condition. Then, M' is comprised 
of edges from left-unmatched vertices to right-unmatched vertices of M, and M U M' is a matching 
of B{A, B); further, since by hypothesis the matching M' has no right-unmatched vertices, neither does 
M U M'. Consequently, by Theorem [T]-(3), this implies that the system is structurally controllable. ■ 

Note that using Lemma we conclude that the system depicted in Figure [^(a) is structurally con¬ 
trollable, yet by using the characterization in O, it is not possible to obtain the same conclusion. 
Further, Lemma provides only a sufficient condition for structural controllability. Nonetheless, these 
conditions can be verified in a distributed manner in the class of interconnected dynamical systems 
formally introduced next. 

Definition 6: We say that an interconnected dynamical system {A, 5) as in ([T]) is a serial system if 
each vertex of the condensed graph V*{A) has at most one outgoing edge. o 

Although serial systems seem a restrictive class of systems, they may exhibit a rich structure as the 
condensed graph of a serial system in Figure [^illustrates. Further, recall that serial systems enable us to 
verify the sufficient conditions for structural controllability in Lemma [^ in a distributed manner. Thus, 
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Fig. 3: In (a), we present the digraph associated to an interconnected dynamical system, where the different subsystems are represented 
inside the dashed boxes. Recall the definition of cactus, it can be readily be seen that the digraph I^(A, B) is spanned by the input cactus, 
depicted by the red edges, rendering the structural system {A,B) structurally controllable by Theorem [^(2). In (b), however, we depict 
possible cacti that span each of the subsystem digraphs. Since a spanning cactus for B{A 2 , B 2 ) has to include a stem comprising at least 
a vertex, neither of the cacti spanning B[A\, B\) and T>{A2, B2) can contain any cycles, and there is no way of prolonging the stem that 
spans T>{A\, B\) to include a stem spanning T>{A2, B2)- This shows that the conditions proposed in |8| are not necessary. 



Fig. 4: Example of a possible condensed graph for a serial system, where each vertex represents a subsystem, and each directed edge a 
non-zero connection matrix, see Definition 

in Algorithm we present the proeedure that eaeh agent employs in order to verify the eonditions in 
Lemma |3l 

Before introdueing Algorithm we explain the funetions that are used throughout the algorithm. All 
these funetions should be able to be applied by the i-th subsystem, Send(t, j) sends the value of the 
variable x to the j-th subsystem, while Rcv(j) makes the system wait to reeeive a message from the 
j-th subsystem and subsequently reads this message. Note that both these funetions ean only be applied 
when systems i and j internet with eaeh other. For eommunieation to be sueeessful, we assume that the 
systems perform these steps synehronously (i.e. they wait for the responses of their neighbors). Finally, the 
proeedure MinWtMaxMatch(c, G) ealeulates the minimum weight maximum matehing on the graph 
G using the eost funetion c, and boolean expressions eontained in square braekets get evaluated (to True 
or False). 

Remark 3: The Algorithm ean be easily adapted to eover the ease where eaeh subsystem only has 
one ineoming neighbor. In this ease, the only differenee is that instead of having Bi as the bipartite graph 
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Algorithm 1 Distributed algorithm to verify sufficient conditions given in Lemma for an arbitrary serial system. 

1 : procedure SEQSxRxCTLf^i, Bi, r, 7 ^ 0 ) 

> r — total number of subsystems, Eij = connection matrices 
2 : nghl(i) ^ {j : Eij 7 ^ 0 } 

3: nghO(i) ■<— {j : Ej^i / 0} 

4: nghbs(i)nghl(i) U nghO(z) 

[> send the dynamic matrix to (the unique) outgoing neighbor 
5: for all j G nghO(i) do 

6 : SEND(Ai,ji) 

7: end for 

> receive the dynamic matrix of the incoming neighbors j G nghl(i) 

8 : for all j G nghl(i) do [> nghl(i) = {ji,. .. ,ji} 

9: Aj ^ Rcv(j) 

10: end for 

Ai Ei^j^ . . . 

0 Aj, ... 0 

11 : 

0 0 ... A,, 

12: BJ, 0 ,...,oy 

13: B{A'i,B'i) 

14: function Ci( 2 /, Xfc,) o define the weight function 

15: if y = Xj, and k,j < rii or k,j > rii then 

return 1 
16: else 

return 2 

17: end if 

18: end function 

19: Mi-s-MlNWxMAXMAXCH(ci, Sj) 

20: Wij(Mi) <r- {xj : Xj right-unmatched w.r.t.Mi and j < rii} 

21 : A/” ■<— {xj : Xj state vertex in non-top linked SCC of T>{Ai, Bi)} 

22: mchd(i) [IJR{Mi) == 0 ] 

23: rchcl(i) ^ [Af == 0] 

[> check if the whole system can be structurally controllable according to whether the conditions are satisfied in the current system 
or not 

24: ctld(i) rchd(i) A mchd(i) 

25: for fc = 1,..., r do 

26: for all j G nghbs(i) do 

SEND(ctld(t),y) 
ctld(j) ^ Rcv(j) 

27: end for 

> reconsider the answer in light of the values from the neighbors current answer 
28: ctld(i) •<— ctld(i) /\ ctld(j) 

jGnghbs(i) 

29: end for 

> return True if the system is structurally controllable and False otherwise 

30: return ctld(i) 

31: end procedure 
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associated with the i-th subsystem and all its incoming neighbors, it would be that of the i-th subsystem 
and all its outgoing neighbors. o 

The next result concerns the correctness and complexity of Algorithm 

Theorem 4: Algorithm is correct, i.e., it verifies the sufficient conditions given in Lemma for serial 
systems. Moreover, it has computational complexity O ^ rnax , with Aj = mi + J2je:XiU{i} where 
rrii and rij are the dimensions of the input and state space for the Ath subsystem, and Xj C {1 ,..., r} is 
the set of subsystems which output is provided to the i-th subsystem. o 

Proof: To prove the correctness of Algorithm we start by proving the claim that a minimum 
weight maximum matching M of Bi w.r.t. the weight-function Cj (defined in step [T^ induces maximum 
matchings on B{Ai), as well as on B{Aj) for any subsystem {Aj,Bj) with nonzero connection matrix 
Eij. Let Mi be the matching resulting from restricting M to the edges of B{Ai), and in order to derive a 
contradiction, assume that is not a maximum matching of B{Ai). As a direct consequence of Berge’s 
theorem (see for example Theorem 1 of [l2^ f the set of right-unmatched vertices of any matching contains 
the right-unmatched vertices of some maximum matching, so let M'^ be a maximum matching such that 
C Further, let Si C M he the set of edges from vertices not in Xi to vertices in Xi, and 

let S'i be those edges in Si that end in vertices from Ur^MA- Now, (M\ (MjUS'j)) UM(US'' is a matching 
of Bi with the same number of edges as M (since it has the same number of right-unmatched vertices) 
and with an overall weight lower than that of M (since by hypothesis S'' C Si), which contradicts the 
fact that M is a minimum weight maximum matching. The same argument works for the matching M, 
of B{Aj) with j A f replacing left-unmatched vertices with right-unmatched vertices. 

Now, since {A, B) is a serial system, there is at most one k A ^ with nonzero matrix A/j *. Therefore, 
we let M' and M" be the maximum matchings of B{Ai) resulting from the maximum matchings of Bi 
and Bk, respectively. Then, by Lemma there exists a maximum matching M that has as left-unmatched 
vertices those of M" and as right-unmatched vertices those of M'. Subsequently, we only need to check 
for each subsystem that there is a minimum weight maximum matching of Bi (w.r.t. the weight function 
Wi) that has no right-unmatched state vertices; in summary, in Algorithm we set up the necessary 
structures until step [T^ 


Now, in step 19 the i-th subsystem computes the maximum matching Mj of Bi, and in step 20 the 
system calculates the associated set of right-unmatched vertices. Next, in step the subsystem calculates 
the set of state vertices in a non-top linked SCC of V{Ai,Bi), and in steps 22 and 23 it verifies the 
existence of right-unmatched state vertices of the i-th subsystem w.r.t. the matching Mi, and the existence 
of in a non-top linked SCC of V^Ai, Bi). Finally, the subsystem decides if the whole system is structurally 
controllable or not in steps 2^29 More precisely, after an initial conditions has been chosen and stored 
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in ctld(7), the subsystem updates this variable with the corresponding variable of its neighbors, and 


repeats this r times. Note that after k iterations of the steps 26-29 the subsystem has updated ctld(7) 
with the corresponding values of all subsystems at k edges of distance from it. Since the condensed graph 
of the systems is weakly connected, and there are only r subsystems, ctld(7) = True if and only if all 


subsystems had initially ctld(j) = True. Finally, in step 30 the subsystem returns the value True or 
False depending on whether or not the system satisfies the conditions of Lemma 

Lastly, the complexity of Algorithmis given as follows: since all of the steps have linear complexity 


except determining the minimum weight maximum matching of Bi in step 19 for which the Hungarian 
algorithm can be used with complexity O (|Ajp), with Ni= pi + ^ {O’ ^ 


{0 ,and Xj C {1 ,... ,r} is the set of indexes of subsystems incoming to the f-th subsystem 
This procedure has to be applied to each of the r subsystems, which implies that the complexity of the 


algorithm becomes O ( max Nf ). ■ 

Remark 4: If the system was not serial then there could be a subsystem, k which outputs to both the i- 
and j-th subsystems. This means that while computing maximum matchings of Bi and Bj separately, we 
can match the state vertex of the /c-th subsystem to two different state vertices, one of the z-th subsystem 
and one of the j-th subsystem. Further, if there is a subsystem with incoming edges from every other 
system, the algorithm will calculate a maximum matching in a centralized manner. 

Now, we move toward distributed algorithms that can verify structural controllability of interconnected 
dynamical systems at large. In this, each subsystem is required to share only partial information about 
its structure with its neighbors. This algorithm, however, has a higher computational complexity than 
Algorithmic In order to infer structural controllability, we apply Theorem [^-(4), and begin by presenting 
an algorithm to verify if each of the state vertices, in the digraph associated to an interconnected dynamical 
system as in ([T]), has a path from an input vertex to it. 

Theorem 5: Algorithm is correct (i.e., it returns True if and only if every state vertex in the z-th 
subsystem digraph is input-reached). Further, Algorithm has complexity 


O ^max Nr, N max ’ 


where n, is the dimension of the state space of the z-th subsystem, and N = X]I=i where fc* is the 
number of SCCs in the z-th subsystem digraph. 


o 


Proof: Note that, since each subsystem can establish two-way communication with its neighbors, 
the communication graph is strongly connected, and thus the instructions in steps [Tj-fTS] only need to be 
executed (at most) r times in order to receive all pairs {id, #SCCs) in the system. Subsequently the total 


number of SCCs, N, can be computed in step 14 
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Algorithm 2 Distributed algorithm to verify condition (4i) of Theorem 

1: procedure REACHED(Ai, Bi, ^ 0, Ek,i / 0, r) 

2 : nghl(i) {j : Eij 7 ^ 0 } 

3: nghO(i) ■<— {j : Ej^t ^ 0} 

4: nghbs(i)-4—nghl(i) U nghO(i) 

5 : ^ #{SCCs of ©(Ai)} 

6 : SCCs(i) 4-{(i, A^i)} 

[> the subsystems interact with each other to learn how many SCCs each subsystem has, in order to find the necessary number of 

communication steps 
7: for fc = 1,..., r do 

8 : for all j G nghbs(i) do 

9: SEND(SCCs(i), j) 

10: SCCsO') ^ Rcv(j) 

11: SCCs(i) 4-SCCs(i) U SCCs(j) 

12: end for 

13: end for 

14: Ar^E;=iSCCs(j) 

15: rchd(i) 4 — {} [> initialized the list of input-reached vertices 

0 add the vertices with incoming edges from input vertices 
16: for j = 1,..., Hi do 

17: if 3fc : {Bi)j^k = 1 then 

18: ADDTo(a;j, rchd(i)) 

19: end if 

20: end for 

21: for k = 1,..., N do 

> transmit to outgoing neighbors which vertices that interact with them have been input reached 

22: for all j G nghO(i) do 

[> M,j is the i-th column of M 

23: SEND({xi : G rchd(i) and {Ej^i),^i 7 ^ 0}, j) 

24: end for 

> add vertices reached from the neighbors’ input reached vertices 

25: for all j G nghl(i) do 

26: avail(j) 4— Rcv(j) 

27: rchd(i) 4— rchd(i) U {xt : {Eij)t,i = l,xi G avail(j)} 

28: end for 

> verify which vertices are reachable from the inputs by using k edges between subsystems steps 

29: for i = 1,..., rii do 

30: rchd(i) 4— rchd(i) U {xt : (Ai)t,s = l,Xa G rchd(i)} 

31: end for 

32: end for 

0 return True if every state vertex the f-th subsystem digraph is input-reached, and False otherwise 
33: return [#rchd(f) == rii] 

34: end procedure 
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Now, assume that each subsystem has a strongly connected state digraph V{Ai). Then, if the system 
has an input vertex, i.e., if Bi A 0, each of the state vertices of the i-th subsystem is added to rchd(i) 
in the first iteration of the for-loop in steps 21 namely in the for-loop [29| - f3T| Further, note that in the 
case where each subsystem has a strongly connected system digraph, N = r, and a path from an input 
vertex to a state vertex contains at most r edges between different subsystem digraphs. Therefore, in this 


case, after N iterations of steps 21-31 all vertices that may be reached by a path from an input vertex 
have been added to rchd(i). 

Alternatively, if the i-th subsystem is not strongly connected, then, assume, without loss of generality, 
that Ai is a block matrix, with submatrices A\^... ,A\ along the diagonal so that V{A\)^ ..., V{A^j) are 
strongly connected. Further, let B},..., B^ be the restriction of Bi to the rows in used by A},, A\, 
respectively. Then, consider the interconnected dynamical comprising comprising, instead of the i-th 
subsystem (A*, Bi), the subsystems [A], B }),..., {A\, Bl) connected amongst them and to other subsys¬ 
tems according to A*. By applying this procedure to every subsystem whose state digraph is not strongly 
connected, we obtain an interconnected dynamical system, where each subsystem has a strongly connected 
digraph. Note also, that we did not change the state digraph of the overall system, thus, a state vertex in the 
overall system digraph is input-reached if and only if it was input-reached in the original system digraph. 
Now, since this system has N = ^{SCCs in all subsystems of the original system digraph} subsystems, 
from the previous paragraph we conclude that after N iterations of steps 21-^ every state vertex in 
P(Ai) that is input-reached in V{A,B) has been added to rchd(A- 

Thus, we have proven that for any interconnected dynamical system, every state vertex of V{Ai) has 
a path from some input vertex in the overall system if and only if #rchd(A = n*. 

Finally, we analyze the complexity of Algorithm]^ the SCCs of T){Ai) can be computed in 0{ni), 
and each of the steps in the for-loop |7-13 can be executed in constant time, which implies that the 


for-loop incurs in complexity (T (r#nghbs(A) which is bounded by C>(r^). Further, the steps 22-24 


and 25-28 can be executed in constant complexity, thus these loops incur in complexity O (#nghO(A) 
and ( 7 ^nghl(A), respectively. Finally, the for-loop in steps [29l - [3n incurrs in linear complexity (on the 
number, Ui, of state variables). So in conclusion, the complexity of Algorithm becomes 


O max < r , Nr, N max Ui 


Next, we present a distributed algorithm to verify structural controllability when the subsystems only 
have access to neighboring subsystems. Briefly, the algorithm verifies both conditions (4?) and [An) of 
Theorem in a distributed manner: on one hand, the condition (4*) of Theorem can be verified by 
applying Algorithm]^ On the other hand. Theorem [^(4u) requires one to compute a maximum matching 
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in a distributed manner. This can be achieved by reducing the problem of finding a maximum matching 
to that of computing a maximum flow [l251l . However, since we only need to detect the existence of right- 
unmatched vertices, we only need to compute a maximum preflow, which corresponds to a flow, where 
the flow on the incoming edges need not be equal to the flow on the outgoing edges of each vertex. So, 
we employ the distributed algorithm provided in [|^ . In order to achieve this reduction, one takes the 
overall system bipartite graph and provides an orientation to each edge, from left-vertex to right-vertex. 
Then, adds two extra vertices, called source and sink. Finally, one adds an edge from the source to each 
of the left-vertices of the bipartite graph, and from each of the right-vertices to the sink and assigns to 
each vertex a capacity of 1 [|^ . The computation of the maximum flow is then done distributedly, where 
each subsystem works to maximize the flow from the source to the sink within a region of the graph 
comprising the subsystems bipartite graph, the source and the sink (note that the source and sink lie in 
all regions, which does not impair the distribution of the algorithm, since the systems need not keep track 
of the excess on the source or the sink), and any vertices in other subsystems to which the system is 
connected. This is achieved through a push-relabel algorithm, briefly described as follows: each of the 
vertices in a region keeps track of an excess (which corresponds to the difference between the incoming 
and outgoing flow), and a label or height. The excess is then pushed from higher labels to lower labels 
increasing the flow through the edges between them until it reaches the sink, or the boundary. Once this 
is achieved, the excess accumulated in the boundary is passed to the corresponding neighboring region, 
and the iterations begin again. However, the existence of boundary vertices limits the parallelization, as 
two instantiations of the algorithm can only (in general) be computed simultaneously, if the regions do 
not share vertices other than the source or the sink. 

From this point onwards we refer to the individual instances of the parallell region discharge algorithm 
presented in [|26l, by PRD. Further, we assume PRD takes as parameters the digraph on which it operates, 
the capacity function, and the neighbors with which it shares vertices other than the source or the sink, 
and returns a maximal preflow on the digraph. 

Theorem 6: Algorithm is correct, i.e., it verifies Theorem [^(4) in a distributed fashion. Further, it 
has a computational complexity of 

O I maxjr^, A^r, max ni,rB‘^ max nf} I , 

V *=1’--’' J 

where B is the number of boundary vertices, and the remaining variables are the same as described in 
Theorem [51 o 

Proof: In order to verify the correctness of Algorithm]^ we have to check if both conditions (4i) and 
{Ail) of Theorem are verified. Furthermore, in order to perform this verification in a distributed manner, 
each subsystem must verify that all vertices in its digraph are input-reached in V{A,B), which is done 
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Algorithm 3 Distributed algorithm to verify condition (4) of Theorem 

1: procedure CONTROLLED(Ai, Bi, Ei^k ^ 0, ^ 0, r) 

2 : nghl(i) {j : Eij ^ 0 } 

3: nghO(i) ■<— {j : Ej^i ^ 0} 

4: nghbs(i) -4— nghl(i) U nghO(i) 

0 verify if every state vertex is input-reached by deploying Algorithm 
5: rchd(i) 4- REACHED(Ai, Bi, Ei^u + 0, Ek,i 7 ^ 0, r) 

[> we set up the graph for applying the Parallel region discharge, where s and t correspond to the source and sink, respectively, the x 
and u vertices correspond to state and input vertices. The upper index i is the index of the subsystem they belong to, and the upper 
index R and L indicate if they are right or left vertices 

7: ^ : {Ai)j,j = l}u{{n],xif) : (BOp,, = 1} 

8 : for all j G nghbs(i) do 


9: 

Vi,j G- 


: (Ej.i) 

h,» / 0} 


10: 

^i,j ^ 


,xi’^} 

• {Ej,i)k,l - 

= 1} 

11: 


{xt^ ■■ 

■■ (E,,r) 

m,k 7^ 0} 


12: 

Sj^i ■<— 

Uxi:^ 

i,R\ 

,x{ ) 

■ {Bi,j)l,k - 

= 1} 

13: 

end for 





14: 

Ss ^ {s} 

xK" 

‘ \ '^i 

Jfc = l 




15: Si ^ X {0 

16: s ^ s^vj St yj Sr,iVj U {Sj,iyjSj,i) 

_jGnghbs(i) 

17: V^ViU U {Vj,iUVr,j) 

jGnghbs(i) 

18: V^{V,S) 

19: function c(e G Si) 

20: return 1 t> all edges have unitary capacity 

21: end function 

0 Deploy a Parallel Region Discharge algorithm to obtain a preflow / on D, with capacity function c 
22: / G-PRD(D, c, nghbs (i ) ) 

23: mchcl(i) 4- [ /(e) == ^i] 

eGSt 

0 check if the whole system can be structurally controllable according to whether the conditions are satisfied in the current system 
or not 

24: ctld(i) 4— rchd(i) A mchd(i) 

25: for fe = 1,..., r do 

0 reconsider the controllability of the overall system, in light of the data received from the neighbors 
26: for all j G nghbs(i) do 

27: SEND(ctld(i), y) 

28: ctld(y) ^ Rcv(j) 

29: ctld(i) 4— ctld(i) A ctld(y) 

30: end for 

31: end for 

32: return ctld(i) 

l> return True if the system is structurally controllable and False otherwise 
33: end procedure 
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by employing Algorithm in step and that none of its state vertiees are right-unmatehed w.r.t. some 
maximum matehing of B{A,B). In addition, reeall that it was already argued in the proof of Theorem]^ 
that the for-loop in steps [25| - f3T] determines if these eonditions are violated in any of the subsystems. 

Now, to verify that Algorithm determines if there are right-unmatehed vertiees in the i-th subsystem, 
we note that in steps [6}-[2T] we generate the digraph V eomprising the right- and left-vertiees of the i-th 
subsystem bipartite graph, and the boundary vertiees of the z-th region aeeording to the preeepts of [[26]l . 
Onee the digraph V is eomputed, we apply PRD to it in step thus obtaining a preflow from souree to 
sink on V whieh is maximum amongst preflows on the whole graph. By the guarantees provided in [|26l . 
together with the equivalenee between the maximum matehing and maximum flow problems, presented 
in [|^ we guarantee that /(^) equal to the number of right-matehed vertiees in a maximum matehing 
of the system bipartite graph, that are state vertiees of the z-th subsystem. So, by eomparing ^ /(e) with 


e&St 


rti in step 23, we are able to infer if there are right-unmatehed vertiees in the z-th subsystem w.r.t. some 
maximum matehing B{A,B). Thus the algorithm returns True if and only if every state vertex of the 
system digraph is input-reaehed, and there are no right-unmatehed vertiees in the system bipartite graph. 

Now, sinee all of the steps of the algorithm have linear eomplexity exeept for step and step 
the eomplexity of Algorithm is given by the maximum of these. Knowing that step has a eomplex¬ 


ity of O ^max 


Nr, N max n* 


(see Theorem 


where N is the number of SCCs on eaeh of 


the subsystems), all that remains to infer is the eom^exity of step 22 This algorithm, as deseribed 
in [12^ . has eomplexity O (n^) [l27l . Il28ll . and the neeessary iterations of the region diseharge that 
eaeh subsystem has to eomplete, ean be bounded by where f3 is the number of boundary vertiees 
in the whole system bipartite graph (that is, the number of vertiees in the bipartite graph that have 
to be shared by several subsystems). Also, in the worst-ease seenario, where eaeh of the subsystems 
is eonneeted to every other subsystem, the region diseharge steps have to be exeeuted sequentially. 
So, the eomplexity of step 22 is given by O ( r/3‘^ max nf ), resulting in an overall eomplexity of 




(9 I max < A^r, max {zzj,r/3^n?} 


2 = 1 ,...,r 


V. Illustrative Examples 

Let us eonsider the eontrol system, whose system digraph is depieted in Figure |^1. This eonsists 
of four intereonneeted subsystems, whose state digraphs are enelosed by grey dashed boxes. We would 
like to assess if this system is strueturally eontrollable from the single eontrol input zzi, eonsidering 
only the loeally available information, i.e., in a distributed fashion. To solve this problem, we apply 
Algorithm (and Algorithm whieh is required as a subroutine) that is exeeuted on subsystem level, 
henee, emphasizing its distributed nature. 
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Fig. 5: Example of the procedure of Algorithm [^applied to the system digraph presented in Figure 
(depicted inside of the dashed boxes), only one of which has an input edge (labeled ui). In each subfigure, the blue edges represent those 
that comprised a path from an input vertex, and the green edges denote those that were added in this iteration or communication step of 
the algorithm. Finally, the even-labeled subfigures correspond to an iteration of Algorithm and the odd-labeled ones correspond to a 
communication step between subsystems. 


5 -1, comprising 4 different subsystems 


In Figure [^1, we present the digraph assoeiated to an intereonneeted dynamieal system eomprising 
four subsystems. Sinee only subsystem [Ai^Bi) has an input vertex, it readily follows from Theorem 
that none of the other subsystems ean be strueturally eontrollable. Now, we deploy Algorithm to verify 
the struetural eontrollability of the intereonneeted dynamieal system. After the initialization steps are 


eompleted, we deploy Algorithm the iterations of whieh can be seen in Figure |^2 to Figure [^8: in 
each iteration (even-labeled subfigures) new vertices are seen to be input-reached (the targets of the green 


edges), and in each communication step (odd-labeled subfigures) the subsystems interact to its outgoing 


neighbors which of their vertices are reached after the iteration has been completed. As can be seen in 
Figure [^8, all vertices have been reached after four iterations, which in this case, corresponds to the 
number of SCCs in all subsystems. 

In Figure we consider an example of a run of the region discharge algorithm running on the bipartite 
graph associated to the digraph in Figure [^1. In this example, we begin applying PRD in step 1 by 
initializing the labels at 2 for each left-vertex, and at 1 for each right-vertex; we also saturate all edges 
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Fig. 6 : Example of the PRD algorithm applied to the verification of structural controllability of the system presented in Figure 5l-l. For 
convenience of referencing, the vertices are given labels rather than colors (X 4 being the black vertex in each subsystem and the others can 
be easily inferred). At the left of each left-vertex and at the right of each right-vertex we insert two numbers, the one in green represents 
the excess of the corresponding vertex, whereas the red number represents its label. Edges in red represent those where the capacity has 
been saturated; and right-vertices in red represent the ones for which the edge to the sink has been saturated. Finally, the vertices in blue 
represent vertices that belong to other regions, i.e., boundary vertices. In order to simplify, in this, we do not include boundary vertices from 
incoming subsystems. Note also, that in this instance, we can run the algorithm in all regions simultaneously, since by not considering the 
incoming edges from other regions in the region graph, we do not allow for flow to be sent back through these edges. 


from the source, which makes it so that all of the left-vertices start with an excess of 1. By successive 
pushing and relabeling, we reach the configuration in 2 where all the excess has either been pushed to 
the sink (and thus the corresponding right-vertex is presented in red) or to the boundary of the region. 
In step 3, we discharge the excess from the boundary into the adjacent regions so that, for example, the 
right-vertex xs in each of the regions has now an excess of 1. Finally, by applying push-relabel again 
in each of the regions, we reach step 4 where all the edges from right-vertices to the sink have been 
saturated (and are thus displayed in red) showing that there is a maximum preflow saturating all edges 


to the sink, and equivalently that there is a maximum matching with no right-unmatched vertices. So, in 
combination with the analysis of Figure we conclude that the interconnected dynamical associated to 
the digraph system in Figure [^1 is structurally controllable. 


VI. Conclusions and Further Research 

In this paper, we have provided several necessary and/or sufficient conditions to verify structural 
controllability for interconnected linear time-invariant dynamical systems based on the local information 
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accessible to each subsystem. Subsequently, we have provided distributed and efficient (i.e., polynomial in 
the dimension of the state and input) algorithms to verify a necessary and sufficient condition for structural 
controllability. The results presented readily extend to discrete time-invariant interconnected dynamical 
systems, since the controllability criterion stays the same. Further, by duality between controllability and 
observability the results also apply to structural observability verification of discrete/continuous linear time- 
invariant interconnected dynamical systems. Whereas the results presented pertain to verifying structural 
conditions, it would be of interest to address design problems; for instance, which state variables need 
to be actuated, or which inputs should be used, to ensure a given structural property. On the other hand, 
it would be of interest to understand if the conditions provided could be adapted to the case where 
some of the entries in the structure of the subsystems and their interconnections are known exactly 
(which corresponds to the case where only some of the components of the overall system are assumed 
to be reliable). Ultimately, such an extension would shed light on the relationship between structural and 
non-structural system-theoretic properties; hence, leading to a better understanding of the resilience and 
performance of interconnected dynamical systems. 

VII. Acknowledgements 

The first author would like to thank ACCESS Linnaeus Center, for their hospitality, as much of the 
writing of this paper and the research that led to it were carried out while visiting KTH in the Spring of 
2014. 


Reeerences 

[1] U. Ozguner and H. Hemani, “Decentralized control of interconnected physical systems,” International Journal of Control, vol. 41, 
no. 6, pp. 1445-1459, 1985. 

[2] E. Davison and U. Ozguner, “Characterizations of decentralized fixed modes for interconnected systems,” Automatica, vol. 19, no. 2, 
pp. 169 - 182, 1983. 

[3] E. Davison, “Connectability and structural controllability of composite systems,” Automatica, vol. 13, no. 2, pp. 109 - 123, 1977. 

[4] J.-M. Dion, C. Commault, and J. V. der Woude, “Generic properties and control of linear structured systems: a survey.” Automatica, 
pp. 1125-1144, 2003. 

[5] C. Lin, “Structural controllability,” IEEE Transactions on Automatic Control, no. 3, pp. 201-208, 1974. 

[6] R. W. Shields and J. B. Pearson, “Structural controllability of multi-input linear systems,” IEEE Transactions Automatic Control, vol. 
AC-2L no. 3, 1976. 

[7] B. D. O. Anderson and H.-M. Hong, “Sbuctural controllability and matrix nets,” International Journal of Control, vol. 35, no. 3, pp. 
397-416, 1982. 

[8] C. Rech and R. Ferret, “About structural controllability of interconnected dynamical systems,” Automatica, vol. 27, no. 5, pp. 877 - 
881, 1991. 

[9] K. Li, Y. Xi, and Z. Zhang, “G-cactus and new results on structural controllability of composite systems,” International Journal of 
Systems Science, vol. 27, no. 12, pp. 1313-1326, 1996. 



26 


[10] L. Blackball and D. J. Hill, “On the structural controllability of networks of linear systems,” 2nd IFAC Workshop on Distributed 
Estimation and Control in Networked Systems, pp. 245-250, 2010. 

[11] G.-H. Yang and S.-Y. Zhang, “Structural properties of large-scale systems possessing similar structures,” Automatica, vol. 31, no. 7, 
pp. 1011 - 1017, 1995. 

[12] S. Pequito, S. Kar, and A. P. Aguiar, “A framework for structural input/output and control configuration selection of large-scale 
systems,” To Appear in IEEE Transactions on Automatic Control, 2015. [Online]. Available: http://ieeexplore.ieee.org/stamp/stamp. 
jsp?tp=&amumber=7112630&tag=l 

[13] S. Pequito, S. Kar, and A. Aguiar, “On the NP-completeness of the constrained minimal structural controllability/observability 
problem,” in To Appear in Automatica, 2014. [Online]. Available: http://arxiv.org/abs/1404.0072 

[14] S. Pequito, A. P. Aguiar, and S. Kar, “Minimum Cost Input/Output Design for Large Scale Linear Structural Systems,” ArXiv e-prints, 
Jul. 2014. 

[15] T. Zhou, “On the controllability and observability of networked dynamic systems,” Automatica, vol. 52, no. 0, pp. 63 - 75, 2015. 
[Online]. Available: http://www.sciencedirect.eom/science/article/pii/S0005109814005020 

[16] C. Chen and C. Desoer, “Controlability and observability of composite systems,” IEEE Transactions on Automatic Control, vol. 12, 
no. 4, pp. 402-409, August 1967. 

[17] W. Wolovich and H. Hwang, “Composite system controllability and observability,” Automatica, vol. 10, no. 2, pp. 209 - 212, 1974. 

[18] Y. Yonemura and M. Ito, “Controllability of composite systems of tandem connection,” IEEE Transactions on Automatic Control, 
vol. 17, no. 5, pp. 722-724, October 1972. 

[19] S. Wang and E. Davison, “On the controllability and observability of composite systems,” IEEE Transactions on Automatic Control, 
vol. 18, no. 1, pp. 74-74, February 1973. 

[20] E. Davison and S. Wang, “New results on the controllability and observability of general composite systems,” IEEE Transactions on 
Automatic Control, vol. 20, no. 1, pp. 123-128, February 1975. 

[21] J. P. Hespanha, Linear Systems Theory. Princeton, New Jersey: Princeton Press, Sep. 2009. 

[22] T. H. Cormen, C. Stein, R. L. Rivest, and C. E. Leiserson, Introduction to Algorithms, 2nd ed. McGraw-Hill Higher Education, 2001. 

[23] J. Munkres, “Algorithms for the assignment and transportation problems,” Journal of the Society of Industrial and Applied Mathematics, 
vol. 5, no. 1, pp. 32 - 38, March 1957. 

[24] C. Berge, “Two theorems in graph theory,” Proceedings of the National Academy of Sciences of the United States of America, vol. 43, 
no. 9, p. 842, 1957. 

[25] R. K. Ahuja, T. L. Magnanti, and J. B. Orlin, Network flows: theory, algorithms, and applications. Prentice Hall. 

[26] A. Shekhovtsov and V. Hlavac, “A distributed mincut/maxflow algorithm combining path augmentation and push-relabel,” International 
Journal of Computer Vision, pp. 315-342, September 2013. 

[27] R. K. Ahuja, M. Kodialam, A. K. Mishra, and J. B. Orlin, “Computational investigations of maximum flow algorithms,” European 
Journal of Operational Research, vol. 97, no. 3, pp. 509-542. 

[28] A. V. Goldberg, “The partial augment-relabel algorithm for the maximum flow problem,” in Algorithms - ESA 2008, ser. Lecture Notes 
in Computer Science, D. Halperin and K. Mehlhom, Eds. Springer Berlin Heidelberg, no. 5193, pp. 466-477. 


